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Resumen 


En el área del estratovolcán La Malinche la fuente de agua es el acuífero. 
El objetivo de este trabajo fue analizar la relación entre población y 
niveles estáticos (NE) en los acuíferos Alto Atoyac y Huamantla. Se 


utilizaron las siguientes variables: número de habitantes (NH), tasas de 
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cambio anual de NE, porcentajes de cambio anual de NE (PCAyg), 
tendencias de NE y NH. Los pozos se agruparon con análisis de 
componentes principales (4CP). Los PCAy¿ se compararon con un diseño 
factorial. La tasa de cambio anual promedio de NE fue igual a -0.159 
m:año?* y 6.7% de PCAyg. Las tasas de cambio de NE fueron 
estadísticamente diferentes entre los acuíferos. La relación entre NE y NH 
que resalta fue igual a -16.5 cm:hab"?. Las tendencias de NH fueron 
mayores que las tendencias de NE. El diseño factorial arrojó que entre las 
temporadas de lluvia y estiaje, los PCAyg no tuvieron diferencias 
significativas, pero entre los tipos de concesión sí los hubo. El ACP 
correlacionó 51 pozos con una componente. En resumen, este estudio 
reveló que los PCAyz son mayores en las zonas agrícolas y urbanas. El NE 
en los pozos de uso industrial se abate cinco veces más que en los de uso 
público. La relación entre población y NE fue más clara en la montaña que 
en los valles. El ACP mostró que los pozos de los alrededores de La 


Malinche se diferencian de los demás del área de estudio. 


Palabras clave: uso de agua, regresión lineal, componentes principales. 


Abstract 


In the area of the La Malinche stratovolcano the source of water is aquifer. 
The objective of this work was to analyze the relationship between 
population and water table (WT) in the Alto Atoyac and Huamantla 
aquifers. The following variables were used: Number of inhabitants (NH), 
WT annual change rates, WT annual percentage changes (APCyr), and 
trends of WT and NH. Wells were grouped with principal components 


analysis (PCA). The APC; were compared with a factorial design. The WT 
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annual change rates and APCyy were -0.159 m:year*? and 6.7%, 
respectively. The rates of change of WT, were statistically different 
between aquifers. The relationship between WT and NH that stands out 
was equal to -16.5 cm:hab*?, NH trends were greater than WT trends. The 
factorial design showed that between the rainy and the low water seasons 
there were no significant differences, but between concession type there 
were. The PCA correlated 51 wells with one component. In summary, this 
study revealed that the APC, are higher in agricultural and urban areas, 
the WT one in the wells of industrial use is reduced five times more than 
in those of public use. The relationship between population and WT was 
clearer in the mountains than in the valleys. The PCA showed that the 


wells around La Malinche differ from others in the study area. 


Keywords: Water use, lineal regression, principal components. 
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Introducción 


El agua es un recurso natural insustituible para los seres vivos. El ser 
humano requiere de agua para uso público, agrícola, comercial o 
industrial. Las fuentes del vital líquido son los cuerpos de agua 


superficiales, las aguas subterráneas, la captación de agua pluvial y agua 
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residual tratada para usos que no requieren agua con calidad 


predeterminada. 


A través de la historia, el agua subterránea ha sido una fuente 
primordial para sostener a la humanidad (Fienen € Arshad, 2016). La 
disponibilidad de agua subterránea está disminuyendo conforme las 
entradas bajan y las salidas de agua del acuífero se incrementan por 
extracción intensiva, especialmente en regiones agrícolas áridas y centros 
urbanos (Lall, Josset, 8 Russo, 2020). En regiones tropicales de África, 
Sudamérica y Centroamérica, la extracción de agua subterránea es 
escasa, porque las aguas superficiales son abundantes y no se requiere 


de irrigación en la agricultura (Margat € van der Gun, 2013). 


La extracción global de agua subterránea es de 650 km*-año”?! 
(Margat 8: van der Gun, 2013); el 70% del agua subterránea se usa para 
la agricultura (Fienen 8: Arshad, 2016), y el 24% en la industria (Valencia, 
Mendoza, Vargas, 8 Domínguez, 2006). A nivel global, el agua 
subterránea aporta en promedio el 20% para consumo humano y en 
algunos países puede ser el 90% (Jakeman et a/., 2016). El 50% del agua 
consumida en las zonas urbanas proviene de fuentes subterráneas 
(Fienen 8: Arshad, 2016) y en las zonas rurales aumenta a un valor de 
66% (Carrard, Foster, € Willetts, 2019). 


En el año 1910, México tenía una disponibilidad de agua de más de 
31 000 m*-hab*-año*; en 1950 disminuyó a 18 000 m3-hab"*.año”!; el 
día de hoy es de aproximadamente 3 600 m3-hab"!.año"* (Palma, 2020). 
En la actualidad (año 2021), se desconoce con una aproximación 
aceptable cuántos mexicanos viven escasez extrema de agua. En 


ciudades del norte del país (p. ej., la ciudad de Monterrey), esta situación 
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se hizo evidente en 2022. Breña y Breña (2007) reportaron que, en 2004, 
35 millones de mexicanos se encontraban en escasez extrema de agua y 


43 millones con disponibilidad baja. 


Anualmente, las precipitaciones pluviales aportan a México unos 1 
449.5 km? de agua. De este volumen, se estima que la evapotranspiración 
representa el 72.1%; el 21.4% escurre por los ríos o arroyos, y el 6.4% 
(92.544 km) se infiltra al subsuelo de forma natural y recarga los 


acuíferos (Conagua, 2018). 


En México, el 39.1% de las fuentes de agua son subterráneas y el 
60.9% son superficiales (ríos, arroyos y lagos) (Conagua, 2018; Conagua, 
2019); este nivel de extracción coloca a México como el sexto país 
consumidor de agua subterránea del mundo (Margat 8 van der Gun, 
2013). En la literatura se reportan valores de extracción de agua 
subterránea para México de 27.4 a 25 km3-año”?* (Margat 8: van der Gun, 
2013; Fienen € Arshad, 2016). De este volumen, de 22 a 64% se usa 
para abastecimiento público; de 33 a 72% para irrigación en la 
agricultura, y de 6 a 24% para la industria (Valencia et a/., 2006; Margat 
8 van der Gun, 2013). Sin embargo, Jiménez (2008) reporta que los 
acuíferos aportan cerca de 75% de agua para consumo de la población, 
más aún, para la Zona Metropolitana del Valle de México, este valor 
aumenta a 80.1% y solo en el 19.9% la fuente es de cuerpos de agua 


superficiales (National Research Council, 1995). 


Los principales focos rojos de la disminución del nivel de agua 
subterránea se observan en las partes áridas y semiáridas del mundo, 
como resultado principalmente de la alta densidad de población, fuerte 


dependencia del agua subterránea, pocas y fuertes precipitaciones que 
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tienen bajas tasas de recarga natural que escurren rápidamente (Fienen 
8 Arshad, 2016). Esto conlleva a que en algunos países los acuíferos 
pueden estar sobrexplotados, como es el caso de México (Jakeman et al., 
2016). 


En México hay 653 acuíferos, 105 de estos (16.1%) fueron 
colocados en la categoría de sobreexplotados a finales del año 2017 
(Conagua, 2018). Los acuíferos Alto Atoyac y Huamantla, que 
corresponden al área de estudio de este trabajo, están fuera de esta 
clasificación. Los datos reportados en la literatura (Conagua, 2015a; 
Conagua, 2020a) indican que la disponibilidad de agua en el acuífero Alto 
Atoyac disminuyó 37.1%, pasando de 0.0467 km en 2013 a 0.0294 km 
en 2020. Por su parte, la disponibilidad de agua en el acuífero Huamantla 
aumentó 12.9%, pasando de 0.0139 km* en 2002 a 0.0157 km? en 2020 
(Conagua, 2015b; Conagua, 2020b). En el escenario del área del volcán 
La Malinche, el agua subterránea es la fuente preponderante para uso 
público (ciudad de Apizaco, zona metropolitana de Puebla y Tlaxcala); 
para la agricultura intensiva de riego (valle de Huamantla), y para la 
industria (parques y corredores industriales), lo cual se debe a que en 
esta área no hay fuentes superficiales importantes de agua. El objetivo 
del presente trabajo fue analizar la relación entre la población de las 
localidades y los niveles estáticos de los pozos clasificados por su tipo de 
concesión (público, agrícola e industrial) en los acuíferos Alto Atoyac y 


Huamantla, pues hasta ahora este tipo de estudio no se ha hecho en esta 


zona. 
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Materiales y métodos 


Área de estudio 


Este estudio se realizó en los acuíferos Huamantla 2903 y Alto Atoyac 
2901 en el estado de Tlaxcala. Dichos acuíferos se encuentran en la 
Región Hidrológica del Río Balsas, el primero localizado en la subcuenca 
cerrada de la laguna de Totolcingo y el segundo en las subcuencas de los 
ríos Zahuapan y Atoyac, los cuales desembocan en el océano Pacífico 
(Conagua, 2020a; Conagua, 2020b) (Figura 1). El estratovolcán La 
Malinche es la zona de recarga de tales acuíferos; ambos son de tipo libre, 
heterogéneos y anisótropos, constituidos en su porción superior por 
sedimentos aluviales y fluviales de granulometría variada. La porción 
inferior del acuífero Huamantla está alojada en una secuencia de rocas 
volcánicas y sedimentarias marinas (calizas), y el acuífero Atoyac está 
alojado en una secuencia de rocas ígneas extrusivas de basaltos, tobas y 
andesitas. Más información de los acuíferos se tiene en la Tabla 1. En el 
acuífero Alto Atoyac, año 2011, se registraron 890 aprovechamientos de 
agua subterránea, de los cuales 600 son pozos activos. En el acuífero 
Huamantla, año 2014, se registraron 324 aprovechamientos, de éstos, 


300 son pozos profundos (Conagua, 2020a; Conagua, 2020b). 
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Estados Unidos YN al 
r” AÑ 
México 4 o. de 
É eS Simbología 
dl 3] MA Presa Allangalepec 
A 7 | —Río Zahuapan 
| Valle Puebla-Tlaxcala 
5 EA Estrato volcán La Malinche 
a EM Acuifero Alto Atoyac 
o 250 500 750 1000km a ] = roll 0 10 > 30 
TN 15090 -110%0' -105*0" -100%0" 95% 000" 5 ; : : 
-98%36* -98%24* -98%12* -9890' -97%48* 
Figura 1. Área de estudio, acuíferos en el volcán La Malinche, estado de 
Tlaxcala. 
Tabla 1. Características de los acuíferos Alto Atoyac y Huamantla. 
Parámetro Alto Atoyac Huamantla 
Temperatura media 
16.2 *C 10a16*C 
anual 
Precipitación media 
878 mm 538 mm 
anual 
Región hidrológica No. 18 Río Balsas No. 18 Río Balsas 
, Alto Balsas-cerrada Libres 
Cuenca Alto Balsas-río Atoyac ] 
oriental 
% : Cinturón neovolcánico. Lagos E a 
Fisiografía ] Planicie neovolcánica 
y volcanes de Anáhuac 
a De 0.6 a 11.6 x 10? m*:s* De 2.7 x 10% a 2.4 x 10? m?s? 
Transmisividad 
De 51.8 a 1 002 m*d* De 23.3 a 2 074 m*d* 
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Parámetro Alto Atoyac Huamantla 


Conductividad 


hidráulica 


De0.5a5m:d*! =--- 


Coeficiente de 
De 0.1a1x10? => 


almacenamiento 


Profundidad del nivel De 5a 160 m La Malinche 
estático De 5 a 25 m sur del acuífero 


De 3 a 140 m 


De norte a sur, con . . 
Direcciones norte y sur, 


Dirección preferencial | alimentaciones de los flancos ] : 
convergiendo al oriente a la 


del flujo subterráneo oriental y occidental hacia el 


E laguna de Totolcingo 
acuífero Valle de Puebla 


Afloramientos de rocas 


volcánicas: riolitas, andesitas, 


Geología superficial 
basaltos, tobas y brechas 


volcánicas 


Porción superior del acuífero: 
sedimentos aluviales y 
fluviales de granulometría 


variada. Porción inferior del 


Geología del subsuelo , : 
acuífero: rocas ígneas 


extrusivas, integradas por 
derrames volcánicos, de 


basaltos, tobas y andesitas 


Cuaternario aluvial permeable 


Cretácico, Terciario, (Qal) y volcánico andesítico 


Estratigrafía : oo Aa, 
Cuaternario (Qvab). Terciario andesítico 


basáltico (Tvab). 


Fuente: Conagua (2015a), Conagua (2015b). 
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Número de habitantes 


Los datos del número total de habitantes (NH) de México, del estado de 
Tlaxcala y de las localidades estudiadas en este trabajo se obtuvieron de 
la página del Instituto Nacional de Geografía Estadística e Informática 
(INEGI) (https://www.inegi.org.mx/, fecha de descarga: 2 de marzo de 
2021). En México, se han hecho 16 censos de población de 1900 a 2020; 
con estos datos se obtuvo la ecuación de regresión de NH, en función del 
tiempo, y las tasas de crecimiento de México, Tlaxcala y las localidades. 
Con las coordenadas geográficas se elaboró en QGIS una capa vectorial 
para poner las localidades sobre los acuíferos Alto Atoyac y Huamantla. 


El QGIS es un sistema de información geográfica de licencia libre. 


Tasa de cambio del nivel estático 


Para evaluar la tasa de cambio del NE con respecto al tiempo (años) se 
utilizaron datos de los años de 1997 a 2017 de 73 pozos concesionados. 
Del acuífero Alto Atoyac fueron 50 pozos: 28 para uso agrícola, 4 para 
uso industrial y 18 para uso público. Del acuífero Huamantla fueron 23: 
16 para uso agrícola, 1 para uso industrial y 6 para uso público. El análisis 
de los datos de la fluctuación del NE es el método más directo y simple 
para estimar el volumen de disminución de agua del acuífero (Fienen 8 
Arshad, 2016). Los pozos analizados forman parte de la red de monitoreo 
de agua subterránea que realiza en las temporadas de lluvia y estiaje la 
Comisión Nacional del Agua (Conagua) Dirección Local Tlaxcala. Mediante 


solicitud al Instituto Nacional de Acceso a la Información (INAD), la 


281 


Tecnología y ciencias del agua, ISSN 2007-2422, 


Open Access bajo la licencia CC BY-NC-SA 4.0 Ñ Ñ o E ¿NAS 
(https://creativecommons.org/licenses/by-nc-sa/4.0/) 15(4), 272-317. DOI: 10.24850/j-tyca-2024-04-07 


a 0) Check for updates 
OPEN ACCESS 


Tecnología y 

CienciaszAgua 
Conagua proporcionó los datos de NE de cada pozo. A estos datos se 
calcularon los parámetros como mínimo, máximo, promedio, coeficiente 
de variación, sesgo y curtosis. Se utilizó la prueba de t para evaluar si los 
valores de NE entre los acuíferos Alto Atoyac y Huamantla son 


estadísticamente diferentes. 


La tasa de cambio de la profundidad (m-+año"*) del NE de cada pozo 
con respecto al tiempo (t) se calculó con la pendiente (mps) de la ecuación 
de regresión lineal, para lo cual se utilizó la serie de datos de 1997 a 
2017. La ecuación de regresión lineal se describe como sigue (McBean 81 
Rovers, 1998; Schuenemeyer €: Drew, 2011): 


NE = a+ Mygt + e; (1) 


donde: 
a = ordenada al origen 
e; - N[0,0?] = término de error aleatorio 


La significancia de los valores de my» se estableció a p < 0.05. Para 
visualizar los gradientes relativos de las my, los valores puntuales de 
cada pozo se representaron en isolíneas, para ello se hizo una 
interpolación espacial con la técnica de la ponderación inversa de la 


distancia en el sistema de información geográfica QGIS. 
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Porcentaje de cambio anual del NE 


Un estimador de cambio del NE es el porcentaje de cambio anual (PCAyg), 
el cual se calculó con la ecuación utilizada por Andersen (2019) en 
econometría. En este método no solo se consideran los valores inicial y 
final, sino todos los valores de la serie datos de manera conjunta en el 


término (t — 1) de la siguiente ecuación: 
PCAng = m,(t— 1) x 100 (2) 


donde: 


m1 = pendiente de la ecuación de regresión lineal estimada con la 


siguiente ecuación 
InNNE = b + m;,t+ e; (3) 


InNE = logaritmo natural de NE 
b = ordenada al origen 


En estadística, el empleo de logaritmo natural reduce el coeficiente 


de variación de los valores (Malamud €: Turcotte, 2013). 
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Comparación de PCAyz 


La comparación simultánea y combinada del PCA y, se realizó mediante un 
diseño factorial. La Conagua reporta mediciones del NE hechas en las 
temporadas de estiaje y de lluvias en pozos concesionados para uso 
público, agrícola e industrial. Considerando la temporada de muestreo y 
tipo de concesión o uso de agua de los pozos, se seleccionó un diseño 
factorial de efectos fijos de dos factores. La temporada de muestreo, 
factor A, tuvo dos niveles: lluvia y estiaje. El tipo de concesión del pozo, 
factor B, tuvo tres niveles: público, agrícola e industrial. El modelo 
factorial de dos factores se representa con la siguiente ecuación 


(Schuenemeyer 8 Drew, 2011): 
PCAngjije = 4 +0; +Pj+ (08); + Cri; (4) 
i=1,..,a, ¡j=1,...,b, k=1,..,m 


donde: 
PCAngijx = porcentaje de cambio anual del NE del pozo ijk 
4 = promedio general del PCAys de todos los pozos 
a, = efecto del nivel ¡-ésimo del factor A 
B; 


af;¡ = interacción entre los niveles i y j de los factores A y B 


efecto del nivel ¡-ésimo del factor B sobre el PCAyg 


€x(ij) = error =N (0,0?) 


284 
o = 2024, Instituto Mexicano de Tecnología del Agua. Tecnología y ciencias del agua, ISSN 2007-2422, 
Orem Acces baja le lleanola o BANCA 90 15(4), 272-317. DOI: 10.24850/j-tyca-2024-04-07 


(https://creativecommons.org/licenses/by-nc-sa/4.0/) 


a 0) Check for updates 
OPEN ACCESS 


Tecnología y 


CienciaszAgua 


Evaluación de la relación entre NE y NH 


Para el análisis de la relación de NE y NH por localidad en los aculferos 
Alto Atoyac y Huamantla se calcularon los valores de NH para los años 
intercensales de 2000 a 2020 con las ecuaciones de regresión descritas 


en el apartado “Número de habitantes”. 


La evaluación de la relación entre NE y NH se realizó con la ecuación 


de regresión lineal: 
NH =c+m,NE + e; (5) 


donde: 
c = una constante 
m, = pendiente de la recta que indica la relación entre el NE de un pozo 
NH = localidad cercana de un pozo 


Esta relación se estableció como válida a un nivel de significancia p 
< 0.05. Dado que los censos y conteos de población de 1995 a 2020 
aportan seis datos, se utilizó el segmento recto de la curva exponencial 
de la población en función del tiempo (Figura 2), para calcular el NH de 
las localidades en los años intercensales de este periodo; de esta manera, 
el número de pares de datos para el cálculo de los coeficientes de 


correlación de NH y NE aumentó a 20. 
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Figura 2. Comportamiento del número de habitantes para el Estado de 


Tlaxcala. 


Tendencias de NE y NH 


No fue plausible asociar la tendencia de NH de una localidad dada con la 
tendencia del NE de los pozos cercanos con los datos crudos de estas 
variables debido a que tienen diferentes unidades e intervalos. Por tal 


razón, se estandarizaron los datos de NH y NE con la siguiente ecuación: 


(dato crudo—promedio) 


Dato estandarizado = —. z (6) 
desviación estándar 
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Posteriormente, se obtuvieron las ecuaciones de regresión lineal de 
los datos estandarizados de NH y NE en función del tiempo t. La 


significancia del modelo fue establecida a un nivel de p < 0.05: 
NH =C + Mint + €; (7) 
NE = a + Mint + €; (8) 


Dado que mz tiende a disminuir y my a aumentar con el tiempo, 
los valores de ms se multiplicaron por -1. Las pendientes se 
transformaron a grados de inclinación (4) mediante las siguientes 


expresiones: 
y 180 
¿yg = atan(Myy) * (9) 


¿yg = atan(-1 * myg) * a (10) 
T 


Para determinar si los valores de %4yy Y %yg tienen diferencias 


estadísticas, se utilizó la prueba de t* a p < 0.05, para dos muestras 


pareadas. 
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Análisis de componentes principales 


La combinación lineal de los niveles estáticos de los pozos se estudió 
mediante el análisis de componentes principales (ACP). Una combinación 
lineal es el vector de los niveles estáticos multiplicados por un escalar, lo 
cual indica que los niveles estáticos están relacionados; de otra manera, 
los pozos que no lo están tienen un comportamiento que no cumple con 
este criterio. El ACP es un procedimiento matemático que transforma un 
sistema de variables correlacionadas a otra de variables no 
correlacionadas para reducir su dimensionalidad y determinar las 


combinaciones lineales (Daultrey, 1976; Schuenemeyer € Drew, 2011). 


Los escalares o valores propios (1;) de los componentes principales 
—que a su vez representan las varianzas de los niveles estáticos 
(variables originales) — se calcularon con la ecuación característica de la 


matriz ([R]: 


RA ¿45] = 0 (11) 


donde: 
fRj = matriz de correlaciones de los NE 
A; = valores propios de los componentes principales 
U) = matriz identidad 


Las cargas o correlaciones de los componentes principales con los 


NE se calcularon con la siguiente ecuación matricial: 
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(1) = 18) (Ap (112) 
donde: 

(Ly = matriz de las cargas o correlaciones entre los componentes 


principales y los NE 
Ey = matriz de vectores propios asociados con cada valor propio 
fA3) = matriz diagonal de 2;; los elementos fuera de la diagonal son cero 


Las operaciones matemáticas se realizaron en Excel, con los 


complementos Solver, MMult y Real Statistics. 


Resultados 


Número de habitantes 


El censo de 2020 reportó que el NH en México fue de 126 014 024 
habitantes y del estado de Tlaxcala 1 342 977. De los años 1900 a 2020, 
la población en México aumentó más de 800% y en el estado de Tlaxcala 
cerca de 700%; en número de habitantes, estos porcentajes son igual a 
112.40 y 1.17 millones de habitantes, respectivamente. Hasta 1940, el 
número de habitantes se mantuvo constante, incluso con un descenso por 
la Revolución Mexicana de 1910 a 1924. El comportamiento de 1900 a 
2020 del número de habitantes (NH) en México y en el estado de Tlaxcala 


se describió con un modelo exponencial NH =C+:e** (p < 0.05), donde C 
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y k son constantes, con un comportamiento recto en los últimos 30 años; 
en este lapso, el NH aumentó 55.1% (44.7 millones) y 76.4% (581 mil), 
en México y en el estado de Tlaxcala, respectivamente. En particular, para 
el estado de Tlaxcala, de 1900 a 1910, la tasa de incremento fue de 1 
285 hab:año"?* y para 2010 a 2020 este valor aumentó a 17 304 hab:año” 
1 (Figura 2). 


Tasa de cambio de NE 


La estadística básica mostró que los pozos del acuífero Alto Atoyac 
tuvieron valores de NE promedio, mínimo y máximo iguales a 57, 0.9 y 
184 m, respectivamente. Para los pozos del acuífero Huamantla, los 
valores respectivos fueron de 68, 2.3 y 153 m. Los coeficientes de 
variación y de sesgo fueron mayores en los pozos del acuífero Alto Atoyac 
(70% y 0.8) que en los pozos del acuífero Huamantla (60% y 0.2). Los 
valores de curtosis resultaron igual a 0.5 para el acuífero Alto Atoyac y - 
0.6 para el acuífero Huamantla. La prueba de t mostró que los valores de 


NE entre los acuíferos fueron estadísticamente diferentes (p < 0.05). 


El análisis mostró que en la época de lluvia seis pozos tuvieron tasas 
de cambio (mps) positivas y 67 pozos tuvieron tasas negativas. El 
promedio de my fue igual a -0.159 m:año*?*; el mínimo de -0.838 (pozo 
Huamantla 39); el máximo de 0.677 m:año* (pozo Tecopilco 7), y 
desviación estándar de 0.194 m-año"?. Los valores más frecuentes de myz 
se encontraron entre 0.0 y -0.4 m:año"? en 87% de los pozos. El valor del 
coeficiente de correlación r de la ecuación de regresión (NE versus 


tiempo) resultó significativo (p < 0.05) en 57 pozos. En cuanto a la época 
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de estiaje, seis pozos tuvieron valores positivos de myg. El valor promedio 
de my» fue igual a -0.153 m:año"*; el mínimo de -0.820 maño"? (pozo 
Huamantla 39); el máximo de 0.759 m:año? (pozo Tecopilco 7), y 
desviación estándar de 0.197 m:año*. Las tasas más frecuentes se 
encontraron entre 0.0 y -0.4 m:año"! en 85% de los pozos. El valor del 
coeficiente de correlación r de la ecuación de regresión resultó 
significativo (p < 0.05) en 63 pozos. La Figura 3 muestra las isolíneas de 
los valores de my». La gráfica de NE de los pozos clasificados de acuerdo 
con el tipo de concesión o de uso versus tiempo (Figura 4) muestra una 


tendencia lineal de cambio descendente en todos los pozos. 
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Figura 4. Niveles estáticos de los pozos clasificados por su uso versus 


tiempo. 


Porcentaje de cambio anual del NE (PCAyz) 


Los promedios de PCAy* fueron iguales a -6.4% (-0.159 m:año"!) y -7.3% 
(-0.153 m:año*1), en la temporadas de lluvia y estiaje, respectivamente. 


El valor de PCAyz más recurrente fue de O a -20% (de O a -0.4 m:año"!) 


293 
0 = 2024, Instituto Mexicano de Tecnología del Agua. Tecnología y ciencias del agua, ISSN 2007-2422, 
anemias Ban le end PINE 0 15(4), 272-317. DOI: 10.24850/j-tyca-2024-04-07 


(https://creativecommons.org/licenses/by-nc-sa/4.0/) 


a 0) Check for updates 
OPEN ACCESS 


Tecnología y 

CienciaszAgua 
en 80% de los pozos. En seis pozos, los valores de PCAyg fueron positivos. 
Los valores extremos de PCAys fueron de -50% (-0.8 m:año"!) en el pozo 
Huamantla 39 y =50% (0.7 m:año"1) en el pozo Tecopilco 7. La 
distribución espacial de los valores de PCAys se encuentran en las isolíneas 
de la Figura 5. Las dos temporadas, lluvia y estiaje, tuvieron isolíneas 


comparables. 
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Comparación de PCAyz 


El análisis de varianza del diseño factorial de efectos fijos de dos factores 
(Tabla 2) para comparar de manera simultánea y combinada los pozos de 
acuerdo con su tipo de concesión y temporada de muestreo mostró que 
los valores de PCAys en las temporadas de lluvias y estiaje no tuvieron 
diferencias significativas, cuyos valores promedio fueron de -6.9 y -7.5%, 
respectivamente (Figura 6a). Mientras que los PCAys de los pozos 
concesionados para uso agrícola, industrial y público tuvieron diferencias 
significativas (p < 0.05). De acuerdo con la Figura 6b, estas diferencias 
se deben a que los PCAyy en los pozos concesionados para uso agrícola e 
industrial son más pronunciados que en los pozos de uso público. El PCAyz 
promedio en los pozos de uso público fue de -2.4% (n = 24); en los pozos 
de uso agrícola fue de -8.9% (n = 44), y en los pozos de uso industrial 
de -10.3% (n = 5). 


Tabla 2. Análisis de varianza del diseño factorial de PCAyz. 


Fuente de variación S.C. g.l. C.M. Fcal. P 
Factor A (temporada) 5.85 1 5.9 0.03 0.87 
Factor B (concesión) 1427.9 2 714.0 3.44 0.03 
Interacción AB 8.20 2 4.1 0.02 0.98 
Residual 29 026.9 140 207.3 
Total 30 468.9 145 


S.C.= suma de cuadrados 


g.l.= grados de libertad 


C.M.= cuadrados medios 
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Figura 6. PCAy en temporada y tipo de uso. 


Relación entre NE y NH 


De los 73 pozos evaluados, 59 se asociaron con el NH de localidades 
aledañas (Figura 7), con los cuales se calcularon los valores de m, y r 
para los años de 1997 a 2017. Los valores de m,, que indican la relación 
entre NE de un pozo y el NH de una localidad, fueron negativos en 55 
asociaciones y en 4 fueron positivos. El valor negativo de m, más 
pronunciado fue igual a -16.5 cm-hab*! (localidad Ocotitla), seguido de un 
valor de -4.9 en la localidad de Mesa Redonda. En 10 localidades, los 


valores de m, se encontraron de -1 a -2 cm:hab"!; en 23 fue de -1 a -9 
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mm:hab”!, y en 20 localidades fue de -0.07 a -0.9 mm-hab“!, Los valores 
positivos m, estuvieron entre 0.2 y 2.6 mm:hab”?, Los valores de r fueron 
significativos (p < 0.05) en 43 asociaciones de NE con la población y en 


16 no lo fueron; es decir, en 73% de los casos la relación fue significativa. 


Habitantes 
o 10000 20000 30000 40000 50000 -60000 70000 


Su e. ee. ur a? 


e Pozos acuífero Alto: Atoyac 
+ Pozos aculfero Huamantla 


Figura 7. Relación entre NE y NH. 
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Tendencias de NE y NH 


El análisis de las tendencias —negativas de NE y positivas de NH— con 
respecto al tiempo empleando datos estandarizados mostró que en 30 
casos los grados de inclinación de las líneas de tendencia del nivel estático 
(4y) fueron mayores que los grados de inclinación de las líneas del 
número de habitantes (44) y en 25 casos el resultado fue inverso. El 
promedio de los valores de %ys fue igual a 5.2% y el promedio de los 
valores de 4yy fue igual a 6.3%. La prueba t* de dos colas mostró que los 
promedios de 4ys Y 42yy no tuvieron diferencias significativas (p > 0.05), 
y que el promedio de Zyy fue ligeramente mayor (4.2%) que yg. Sin 
embargo, los valores de %y¿ tuvieron cuatro veces mayor variabilidad que 


los valores de Xyy- 


Análisis de componentes principales 


La Figura 8, Figura 9 y Figura 10 muestran el resultado del análisis de 
componentes principales. La distribución geográfica de los pozos 
correlacionados con los componentes a valores de r > 0.8 y r < -0.8 se 
encuentra en la Figura 8. La Figura 9 muestra los primeros 10 
componentes principales con sus respectivos valores de A; que explican la 
varianza de los valores originales y la varianza acumulada en el eje 
derecho. El componente 1 tuvo un valor de A, igual a 42.0, que representó 
57.6% del total de los pozos; los componentes 2, 3 y 4 representaron el 
6.4, 5.8 y 4.9%, respectivamente; los componentes 5 al 20 representaron 


de 2 a 1%. El 91% de los datos originales fue representado con 10 
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componentes y el 100% con 20. Sin embargo, solo tres componentes 
tuvieron cargas o correlaciones de r > 0.8 y r < -0.8 con las variables 
originales (niveles estáticos), y cinco componentes tuvieron correlaciones 
der>0.6yr<-0.6. La Figura 10 muestra las cargas de los componentes 
1 y 2: la Figura 10a para los pozos calificados con base en el acuífero y la 
Figura 10b con base en el tipo de uso. En ambas gráficas se observa una 
concentración de 47 pozos (64.4%) del lado derecho, los cuales tuvieron 
correlaciones positivas (r > 0.6) con el componente 1; en el lado izquierdo 
se encuentran cuatro pozos (5.5%), con correlaciones negativas (r < - 
0.6) con esta misma componente; es decir, con la componente 1 
estuvieron correlacionados 51 (69.9%) de los 73 pozos. Dos pozos (2.7%) 
tuvieron correlaciones (r > 0.6) con la componente 2. Con la componente 
3 estuvieron correlacionados negativamente (r < -0.6) cuatro pozos 
(5.5%). Con la componente 4 estuvieron correlacionados positiva y 
negativamente dos pozos (2.7%). Con la componente 5 estuvo 
correlacionado negativamente un pozo (1.37%), y 13 pozos (17.8%) no 


estuvieron correlacionados con ningún componente. 
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Figura 10. Correlaciones de los componentes 1 y 2. Pozos clasificados 


de acuerdo con el acuifero(a) y el tipo de uso (b). 


Discusión 


Número de habitantes 


El NH de 1900 a 2020 en el estado de Tlaxcala se multiplicó por siete, con 
una tasa de crecimiento anual igual a 1.96%, la cual es ligeramente 
mayor a la tasa de 1.19% para México, que se encuentra casi en la parte 
media del ranking mundial de tasas de crecimiento de 116 países 


reportados por el WorldData.info (https://www.worlddata.info). En esta 
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misma base de datos se reporta que de 2012 a 2021 el NH en México se 
incrementó 11.07%, mientras que en el estado de Tlaxcala el NH de 2010 
a 2020 aumentó 14.8%. El comportamiento del NH en el área de estudio 
se puede describir con un modelo exponencial, como se puede observar 
en la Figura 2. En los últimos 30 años, el crecimiento del NH fue recto, 
sumándose 581 mil habitantes; este comportamiento del NH es 
ligeramente diferente al de México, el cual tiende a estabilizarse, como se 
muestra en la página del INEGI (https://www.inegi.org.mx). Si el NH 
continúa aumentando de forma exponencial en los próximos 20 años, se 
sumarán unos 400 mil habitantes, lo cual traerá consigo mayor demanda 


y presión sobre el recurso hídrico. 


Tasa de cambio de NE 


Las profundidades promedio de NE indicaron que los acuíferos Alto Atoyac 
y Huamantla son comparables con los reportados para el Valle de 
Atemajac y Toluquilla, Jalisco, México (Hernández-Antonio et a/., 2015). 
Los valores mínimos se ubicaron en pozos del Valle Puebla-Tlaxcala y los 
valores máximos en el área del Volcán La Malinche, lo cual hace evidente 
la influencia de la altitud y topografía. Los valores de sesgo y curtosis 
indicaron que los NE tienen distribución normal con base en el criterio que 
establecen McBean y Rovers (1998) para este estadístico, aunque los 
datos están ligeramente sesgados a la derecha (€, < 1) y concentrados 


alrededor de la media (C, = 0). 


El número de pozos con signos positivos y negativos de my» fueron 


iguales tanto en la temporada de lluvias como de estiaje. El promedio, 
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mínimo, máximo, desviación estándar y frecuencia de los valores de myz 
de las temporada de lluvias prácticamente fueron iguales a las de la 
temporada de sequía, lo cual podría ser un indicio de que la entrada de 
agua por lluvias a los acuíferos Alto Atoyac y Huamantla no compensa la 
extracción. Las lluvias en el acuífero Alto Atoyac son de 878 mm anuales 
y en el acuífero Huamantla de 538 mm. El valor negativo más extremo 
de myz fue igual a -0.838 maño? en el pozo Huamantla 39, para riego 
en el Valle de Huamantla, el cual se aproxima al valor de -1 m-año"! para 
la Zona Metropolitana de la Ciudad de México (Carrera-Hernández 8 
Gaskin, 2007), y es comparable al valor reportado por Rahman, 
Kamruzzaman, Jahan, Mazumder y Hossain (2016) para el noroeste de 
Bangladesh. Los valores intermedios de myz (-0.15 maño!) indican que 
el agotamiento del agua subterránea en el área de estudio se encuentra 
dentro de los intervalos reportados en la literatura (Herbert € Dóll, 2019; 
Hu et al., 2019; Joshi et a/., 2021). En particular, para el acuífero Alto 
Atoyac, la Conagua (2015a) reportó una tasa de abatimiento de 0.1-0.2 
m:año"!, en el periodo 1997-2011, y para el acuífero Huamantla una tasa 
de abatimiento de 0.3-0.5 m:año? (Conagua, 2015b). La prueba de t 
mostró que estas tasas de abatimiento son estadísticamente diferentes 
(p < 0.05). En cinco pozos, los valores de my» tuvieron signos positivos; 
estos pozos se encontraron cerca de la presa de Atlangatepec y de la 
laguna de Acuitlapilco, lo cual significa que dichos cuerpos de agua elevan 
el NE en alrededor 0.65 m-año”?. Joshi et al. (2021) reportaron valores de 
elevación del nivel de agua subterránea de 1.28 m-año*?. En la literatura 
se reporta que la elevación del nivel estático está correlacionada con 


corrientes superficiales y precipitación (Weider 8, Boutt, 2010; Dudley € 
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Hodgkins, 2013). En nuestra área de estudio, la cercanía de los pozos con 


los cuerpos de agua puede ser la causa de la elevación del nivel estático. 


Porcentaje de cambio anual del NE (PCAyg) 


Los valores de PCAyg durante el periodo de 1997 a 2017 se encontraron 
en alrededor de 7%. En las zonas agrícolas de los valles de Huamantla y 
sur de la cuenca del Alto Atoyac, así como en la zona urbana de Apizaco, 
los valores de PCAyg fueron de 10 a 35%; estos valores son altos, 
comparados con los reportados en la literatura (0.05-0.4%); valores de 
PCAyg Mayores a 1% se consideran significativos (Taranaki Regional 
Council, 2016). El análisis de varianza de tipo factorial del PCAyy arrojó 
que no existen diferencias significativas entre las épocas del año; este 
resultado confirma que los cambios del NE no están asociados con el 
clima, sobre todo las precipitaciones pluviales (Rahman et a/., 2016; 
Fienen € Arshad, 2016; Kumar, Chandniha, Lohani, Krishan, 8 Nema, 
2018; Lall et a/., 2020). El análisis de varianza mostró claramente el 
efecto del tipo de concesión sobre el cambio de NE; los pozos 
concesionados para uso industrial tuvieron valores de PCAy¿ más 
pronunciados, seguidos por los de uso agrícola y público. Este resultado 
evidencia que el NE en los pozos de las industrias se abate unas cinco 


veces más que en los pozos de uso público. 
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Relación entre NE y NH 


El análisis de la relación (m,) entre NE y NH mostró mayores valores en 
dos pequeñas localidades aledañas: Ocotitla y Mesa Redonda. Una 
posibilidad de que en estas localidades se registraron los mayores valores 
de m, es que están localizadas en la montaña, en la parte alta de la 
cuenca, al norte del acuífero Huamantla, donde posiblemente las entradas 
de agua, como las precipitaciones pluviales, no compensan la extracción 
de agua del acuífero. Por otro lado, los menores valores de m, se 
calcularon en los pozos que se encuentra en los valles de las cuencas del 
río Zahuapan y del valle de Huamantla. En esta área no fue plausible 
evidenciar con claridad el efecto de NH sobre el abatimiento de NE, 
aunque la densidad de población es mayor comparado con la población 
que habita en la parte alta de la cuenca; esto se debe tal vez a que el 
volumen de agua extraída por los pozos es compensado por las entradas 
de agua superficiales y subsuperficiales; ello se argumenta con el 
resultado obtenido en el pozo Xicohtzingo, que se encuentra en la 
confluencia de los ríos Zahuapan y Atoyac, en el cual se calculó un valor 


de m, cercano a O y un valor de r =-0.06 (p = 0.7922). 


Tendencias de NE y NH 


Las tendencias de NE y NH medidas con los ángulos de inclinación (4) de 
las rectas de regresión mostraron que en 55% de los casos los ángulos 
de NE fueron mayores que los ángulos de NH y en 45% el resultado fue 


inverso. Aunque el ángulo promedio de NH fue ligeramente mayor al 
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ángulo promedio de NE, la prueba de t* de dos colas mostró que los 
ángulos de NE y de NH no tuvieron diferencias significativas (p > 0.05). 
En la literatura se reporta —para las zonas áridas y semiáridas— una clara 
relación del abatimiento del acuífero por mayor demanda de agua debido 
al aumento de la población (Margat 8: van der Gun, 2013; Jakeman et al., 
2016; Islam 8 Islam, 2017; Kumar et al., 2018; Boretti 8 Rosa, 2019; 
Elizondo 8. Mendoza-Espinosa, 2020; Hu et al., 2019). Para el área de 
estudio, la relación del abatimiento del acuífero con el crecimiento de la 
población se observó que no es generalizado en todas las poblaciones y 
pozos, esto se puede deber a que el uso principal del agua es el agrícola 
y a que el área de estudio se encuentra en el altiplano mexicano, donde 
el factor de la altitud y la topografía tienen un rol en la recarga del acuífero 
como se ha reportado en la literatura (Jakeman et a/., 2016; Taranaki 
Regional Council, 2016; Islam 8 Islam, 2017). Este resultado implica que 
hay áreas en la zona de estudio donde hay una relación de NH con NE 
que debe ser tratada de manera diferente de aquellas áreas donde no 


existe, a fin de evitar posible desabasto de agua a la población. 


Análisis de componentes principales 


El análisis de componentes principales mostró que los valores de NE en 
60 pozos son una combinación lineal. Un solo componente representó 
cerca del 60% de la varianza, lo cual indica que NE de los pozos en torno 
a esta componente están en una misma dimensión y su comportamiento 
es colineal. La representación de la varianza aumentó a 75% con cuatro 


componentes. Con la componente 1 estuvieron correlacionados positiva y 
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negativamente 39 pozos (r > 0.8 y r < -0.8), de los cuales 23 se 
encuentran en el acuífero Atoyac y 16 en el acuífero Huamantla; a su vez, 
estos pozos fueron 23 de uso agrícola, 2 de uso industrial y 14 de uso 
público. Los pozos con valores de r < -0.8 asociados con la componente 
1 tuvieron un valor de my positivo. Con la componente 2 se correlacionó 
el pozo Xaltocan del acuífero Alto Atoyac para uso público, el cual se ubica 


cerca del río Zahuapan (Figura 8). 


Es notorio el resultado que de los 33 pozos que no estuvieron 
correlacionados con ningún componente (valores de cero en la Figura 8), 
27 se encontraron en el acuífero Atoyac y 6 en el acuífero Huamantla; 21 
están dedicados a la agricultura, 3 a la industria y 9 para uso público. El 
efecto de la ubicación geográfica de los pozos sobre las correlaciones de 
NE con la componente 1 es evidente debido a que los pozos se 
encontraron en los alrededores del volcán La Malinche y en el Valle de 
Huamantla. Al norte y sur de la cuenca del Atoyac están los pozos que no 
se correlacionaron con ningún componente. Ello sugiere que el análisis de 
la correlación de los niveles estáticos con los componentes principales 
podría servir de apoyo para estudiar los pozos de agua de manera 
holística y mejorar su gestión, a fin de asegurar la provisión de agua a las 


futuras generaciones. 


309 


2024, Instituto Mexicano de Tecnología del Agua. Tecnología y ciencias del agua, ISSN 2007-2422, 


Open Access bajo la licencia CC BY-NC-SA 4.0 Ñ Ñ o E ENA 
(https://creativecommons.org/licenses/by-nc-sa/4.0/) 15(4), 272-317. DOI: 10.24850/j-tyca-2024-04-07 


a 0) Check for updates 
OPEN ACCESS 


Tecnología y 


CienciaszAgua 


Conclusiones 


El comportamiento del número de habitantes de México y Tlaxcala se 
describió con un modelo exponencial, pero en los últimos 30 años fue 
recto. Con este tipo de crecimiento se estima que en los próximos 20 años 
se sumarán 400 mil habitantes en el estado de Tlaxcala, los cuales 


aumentarán la demanda y presión sobre el agua. 


La tasa promedio de NE respecto al tiempo fue igual a -0.159 m:año” 
1. En seis pozos, las tasas fueron positivas y en 67 pozos fueron negativas. 
Las tasas de NE entre los acuíferos Alto Atoyac y Huamantla fueron 
estadísticamente diferentes (p < 0.05). La tasa promedio de NE —en 
términos de porcentaje de cambio anual— resultó igual a -6.4%. Estos 
valores indican que los descensos del nivel estático son comparables con 
los reportados en la literatura para otros lugares y tiempo, pero en 
términos de porcentaje de cambio son mayores. En las zonas agrícolas y 
en la zona urbana de Apizaco, los porcentajes alcanzaron valores de 10 a 
35%. 


El diseño factorial mostró que NE medido en las temporadas de 
lluvias y estiaje no tuvieron diferencias estadísticas, mientras que con 
base en el uso, las diferencias fueron significativas (p < 0.05). 
Gráficamente se observó que los pozos para uso agrícola e industrial 
tuvieron mayores porcentajes de cambio que los pozos de uso público. El 
NE en los pozos de uso industrial se abate hasta cinco veces más que en 


los pozos para uso público. 


La relación entre el nivel estático y el número de habitantes resultó 


significativo en 73% de las localidades estudiadas; el valor más 
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pronunciado resultó igual a -16.5 cm:hab"?! en un pozo y población 
ubicados en la montaña; en 52 casos, los valores se encontraron entre - 
0.01 y -2 cm:hab* en pozos ubicados en los valles. Los valores 
significativos de r se obtuvieron en pozos ubicados en la montaña y los 


no significativos en los valles. 


El promedio de los ángulos de las tendencias de NE fue menor que 
el promedio de los ángulos de las tendencias de NH. Es decir, los valores 
de NE disminuyen con menor rapidez que el aumento de NH. Sin 


embargo, no se encontraron diferencias significativas entre ambos. 


Los niveles estáticos del 57.6% de los pozos fueron representados 
en un componente principal. Los primeros cuatro componentes 
representaron 74.7% de los pozos. No obstante, con tres componentes, 
las correlaciones tuvieron valores r > 0.8 y r < -0.8, y con cinco 
componentes a valores r > 0.6 y r < -0.6. Los pozos correlacionados con 
la componente 1 se encuentran en el valle de Huamantla y alrededores 
de la Malinche. Los pozos que no tuvieron correlación con ningún 
componente están mayoritariamente al norte y sur del acuífero del 


Atoyac. 
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